Quench induced Mott insulator to superfluid quantum phase transition 
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Mott insulator to superfluid quenches have been used by recent experiments to generate exotic superfluid 
phases. While the final Hamiltonian following the sudden quench is that of a superfluid, it is not apriori clear 
how close the final state of the system approaches the ground state of the superfluid Hamiltonian. To understand 
the nature of the final state we calculate the temporal evolution of the momentum distribution following a Mott 
insulator to superfluid quench. Using the numerical infinite time-evolving block decimation approach and the 
analytical rotor model approximation we establish that the one and two dimensional Mott insulators following 
the quench equilibriate to thermal states with spatially short-ranged coherence peaks in the final momentum 
distribution and therefore are not strict superfluids. However, in three dimensions we find a divergence in the 
momentum distribution indicating the emergence of true superfluid order. 
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Introduction Systems of ultra-cold atoms provide us with 
many-body systems with tunable Hamiltonians that may po- 
tentially be used to realize various quantum phases of mat- 
ter. However many of the interesting phases that have been 
predicted require extremely low temperatures that have not 
been reached in experiments. One new phase of matter that 
has been realized in cold atomic systems is the Mott-insulator 
(MI) phase, which is a gapped collective phase and therefore 
has the special property of a low entropy-density. This low 
entropy density makes the MI a favorable starting point for 
creating other low-temperature phases of atoms. In particular, 
a recent set of experiments [ 1 ] uses an effective Mott phase as 
a starting state and suddenly ramps up the tunneling to create 
potentially exotic low temperature superfluid-like (SF) states 
with p and /-orbital symmetries that have never been realized 
before. While such a sudden change or quench of the Hamil- 
tonian from an initially known state to a desired Hamiltonian 
is a promising way to realize new states of matter [2, 3], it is 
not clear whether a system following the quench evolves to a 
state near the desired ground state phase in any sense. Since 
the quenching process is manifestly non-adiabatic and very far 
from equilibrium, it is certainly possible, perhaps even likely, 
that the final state has only a small overlap with the ground 
state of the final Hamiltonian! Measurements of the momen- 
tum density using time-of-flight imaging do indeed show that 
the initially uniform momentum distribution, which charac- 
terizes the MI phase, develops coherence peaks that are char- 
acteristic of the superfluid phase. However, this evolution 
where the momentum density bunches up into peaks instead 
of spreading out is counter-intuitive and suggests dissipation 
of energy similar to evaporative cooling that is used to create 
Bose-Einstein condensates (BECs). Therefore, to understand 
the utility of this approach to realizing new phases of ultra- 
cold atoms, it is necessary to understand both the origin and 
extent of the superfluid state that is formed following similar 
quenches. 

A natural approach to understanding these quenches is to 
consider an analogous quench in the most studied and simplest 
transition that has been observed in ultra-cold atomic systems. 



namely the MI to SF phase transition in s-wave optical lattices 
In addtion there have been several theoretical studies 
of quenches across the MI-SF transition 1.5 - 10,1 which have 
focussed on aspects such as the short-time coherent dynamics 
or the development of local phase coherence and entangle- 
ment in real space which is a particularly interesting aspect of 
MI to SF quenches. However, these studies do not completely 
characterize the final state of the system reached. In particular 
they contain only indirect information about the long-range 
coherence that is measured by time-of-flight imaging. There- 
fore some very recent works lITll [l2t have concentrated on 
comparing the time-of-flight measurements in quench experi- 
ments with theoretical estimates of the momentum density. 

In this letter we develop an understanding of the peak for- 
mation in momentum distribution following an MI-SF quench 
by calculating this distribution for bosonic atoms trapped in an 
s-wave optical lattice. Since the initial MI phase is completely 
incoherent between sites, the momentum distribution is uni- 
form. To reach a SF state, the tunneling is instantaneously 
quenched to a finite value on the SF side of the ground-state 
phase diagram. We find that the final momentum distribu- 
tion after the quench depends strongly on the dimensional- 
ity of the lattice under consideration. The simplest case is 
the zero-dimensional (OD) lattice, where the momentum den- 
sity following the MI-SF quench continues to oscillate in time 
and thus does not appear to equilibrate. In contrast, the one- 
dimensional (ID) MI chain, whose dynamics we calculate nu- 
merically using the systematically exact infinite time-evolving 
block decimation (iTEBD) iflsll approach yields a momen- 
tum distribution which approaches equilibrium. Furthermore, 
the dynamics of the momentum density qualitatively agrees 
with the distribution calculated from the quantum rotor model 
(exact in the large filling limit) |14, 15] studied within the 
semi-classical truncated Wigner approximation (TWA)[16'|. 
Moreover, the equilibrium momentum density estimated from 
iTEBD quantitatively agrees with the result from the classi- 
cal rotor model. Using this fact, which we verified in OD and 
one dimensions (ID), we estimate the equilibrium momen- 
tum density using the classical rotor approximation in two di- 
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FIG. 1 : (a)Quantum dynamics simulations of visibility for OD system 
show absence of equilibration, (b) Equilibrium visibility of momen- 
tum distribution (defined in Eq.[3} following a quench in a OD (2 site) 
Mott insulator from Ji„itiai = to J final = J at various fillings n 
as a function of J /U . Equilibrium is assumed to be thermal. Results 
at large filling n agree well with the rotor value. 



mensions (2D) and three dimensions (3D) following a quench, 
which are relevant to the recent experiments fP] . We find that 
such an MI-SF quench in the 2D limit leads to a classical ther- 
mal gas state with short-ranged momentum correlations (and 
no condensate fraction). On the other hand, an MI-SF quench 
in the 3D case, which directly applies to the experimental re- 
sults, yields a true condensate. 

Quenches in the Base Hubbard model: Ultra-cold bosonic 
atoms in a deep optical lattice can be described by the Bose- 
Hubbard model which is written as 
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where (i, j) are nearest neighbor sites on the optical lattice, fej 

is the boson creation operator on site j, Uj = b^-bj is the num- 
ber operator of bosons on site j and U is the charging energy 
per site. The inter-site tunneling matrix element Q between the 
sites is taken to vanish in the initial state so that the system is 
initially in a MI ground state with {rij) = n (integer) atoms 
per site. The momentum distribution of the bosons, that is 
measured by time-of-flight images, is the fourier transform of 
the off-diagonal density matrix Uk — e**^-' (foj^o). When 
the state following the quench reaches a thermal equilibrium. 



the correlations can be calculated from the free energy 

i^[/3,M] = -ilogrr(e-^^-''^) 



(2) 



where the inverse final temperature /3 and chemical poten- 
tial /i are determined by the conserved initial energy {E) — 
-dpifiF) = and {N) = ~l3df,F = n. 

OD quenches: Let us start by considering the simplest MI- 
SF quench of a OD MI which has two sites connected by a 
bond forming a two-well Josephson junction lHsIl . The Hamil- 
tonian H for such a system can be numerically diagonalized 
exactly to obtain both the equilibrium value and dynamics of 
the momentum distributions Uk, which we plot in Fig.[T] For 
convenience we focus on the "visibility" of the momentum 
distribution. 



Visibility = 



2{n{k = 0) - n{k = ir/a)) 
(n(fc = 0) + n(fc = 7r/a)) ' 



(3) 



which is a measure of how peaked the momentum distribu- 
tion is at its maximum at A: = compared to its minimum 
at fc = ir/a (a being the lattice constant of the optical lat- 
tice). The visibility following a quench in the OD MI-SF 
quench (Fig. ^a)) oscillates as has been found in previous 
results |6, 7] and does not equilibrate in contrast to higher di- 
mensional Mis where we find that reaches the equilibrium 
distribution calculated from Eq. |2] Therefore it is useful to 
consider the results of the thermal distribution in Fig. [Tib). 

Consistent with previous theoretical studies flllslj the re- 
sults in Fig. [T] show that for reasonably large filling factors 
n (other than n = 1), the time-scale for dynamics and max- 
imum visibility attain universal values for intermediate J/U 
corresponding to the quantum rotor model: 
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where TTj = —iQ^^is the angular momentum canonically con- 
jugate to 0j lfT4i [isll and J = 2n(^ is the Josephson coupling. 
The free-energy F is approximated by replacing H — > Hrot 
in Eq. |2l At values of J comparable to the Mott-gap [/, the 
system equilibrates to a sufficiently high temperature so 
that the free-energy F in the classical limit |7] is a combi- 
nation {F = -l2sZxY. 1^) Qf charge fluctuations with 

energy and phase fluctuations described the classical XY 
model 



HxY = - cos ~ 

(ij) 



(5) 



where Zxy = j Yij d(t>je f^^^^ is the partition function of 
the XY model. In this approximation, /3^^ is determined by 
the energy conservation equation 
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where J = /? J and d is the average coordination number of 
each site of the lattice under consideration. For the OD lattice, 
d = \ (after dropping a global U{1) degree of freedom). The 
equilibrium momentum distribution is given by the Fourier 
transform of the real-space correlation function 
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Within the XF-model approximation, the coherent momen- 
tum distribution peak results from a dissipation of phase fluc- 
tuation energy Hxy in the initial state into charge fluctuation 
energy in the final state. The OD MI equilibrium is given by 
J = 1.065 so that one has an equilibrium visibility of — 0.939 
which is in good agreement with the large n equilibrium in 
Fig. [TJb). At large values of J, fluctuations in the number 
density {\{nj — n) — TTj\ > n) invalidate the simple rotor 
approximation iJ^ot ■ 

ID Mott insulators: Unlike the OD MI, the visibility peak 
of the ID MI (Fig.|2|b)), calculated using the in principle ex- 
act iTEBD method, appears to tend to a steady state value. 
These results agree qualitatively with semiclassical dynam- 
ics associated with the quantum rotor model Hrot calculated 
within the TWA I?]. Within the TWA, the distribution of 
the phases 0j begins with the uniform distribution for the MI 
phase and evolves according to the classical equations of mo- 
tion |0,[ll 
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Thus, as seen in Fig.|2jb), the time-scale for dynamics follow- 
ing the MI-SF quench within iTEBD and also the classical 
rotor approximation scales as r oc {JUy^l'^. However, as m 
the OD case, the classical estimate for the final equilibration 
visibility appears to be valid only for intermediate values of 
J/U. 

The partition function for the classical OD XY model 
(Eq. |5]l, can be written in a transfer-matrix form (see for ex- 
ample HI) as ZxY = Tr[M^{<P - 0')] where M((/)) = 
gjcos0 ^ and I„{x) are modified Bessel 

functions. The equilibrium momentum distribution of the 
bosons for the N = (x lattice is calculated to be 



rik 



4J2 + 1 - 4Jcosfc 



(9) 



The classical equilibrium visibility from Eq. |3] at the equi- 
librium parameter J = 1.065, consistent with Fig. |2ja), is 
found to be 1.52. Moreover, the final momentum distribution 
dX n = 1, ^ = Q.bU shown in Fig. |2b) is in quantitative 
agreement with the classical equilibrium Uk- 

2D Mott insulators: The classical rotor model approxima- 
tion can also be used to calculate the equilibrium momentum 
distribution for a MI-SF quench in a 2D MI at filling factors 
71 3> 1 in the initial state. Such a 2D MI quench could be 
realized in experiments by adding a planar confining poten- 
tial to the set-up in recent experiments [1]. Here we use the 




FIG. 2: (a) Quantum dynamics (iTEBD) simulations of visibility for 
ID system and TWA show equilibration following a quench in a MI 
from Jinitiai = to J final — J 3.1 various fillings n and J/U. 
The thermal equilibrium visibility determined from the XY model 
is in good agreement with the iTEBD and TWA results, (b) Time- 
evolution of the momentum distribution for J — 0.5 U and n = 1 
shows evolution of the momentum distribution to one that is in good 
agreement with the classical thermal distribution. 



cluster- update Monte-Carlo (MC) scheme 111811 to calculate the 
average energy for the 2D XY model (Eq.|5} on a 127 x 127 
lattice with periodic boundary conditions to determine the 
equilibrium final temperature satisfying Eq. |6] We find the 
equilibrium temperature = 1.51Jtobe larger than the 



Kosterlitz-Thouless temperature P^^Tp = 0.89 J ||I9|] . As seen 
in Fig. [3] the momentum density at J = f3J = 0.662 cal- 
culated using MC shows a peak (Fig. O at fc = which is 3 
times higher than the ID case and is qualitatively similar to the 
momentum distributions observed in experiment for p— wave 
superfluids HI. However, the equilibrium momentum distri- 
butions of the quenched MI both in ID and 2D remain finite 
and smooth around fc ^ indicating a boson correlator that is 
exponentially decaying in real space. 

3D Mott Insulator: The classical rotor approximation can 
also be applied to the 3D case as well. In contrast to the 
lower dimensional MI-SF quenches, we find that in 3D, the 
energy density e„<; = —0.988 at MI-SF transition calcu- 
lated using MC [23] is slightly higher than the equilibrium 
value following the quench given by Eq. |6l indicating a equi- 
librium phase with a true SF component. This is consis- 
tent with the high-temperature expansion for a 3D MI on a 
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FIG. 3: (a) Equilibrium momentum density as a function of momen- 
tum kx after an MI-SF quench as compared with the initial distribu- 
tion for OD, 2D and 3D Mis. All calculated within the classical rotor 
model approximation. The 3D MI case shows the emergence of a 
true condensate with a divergent fc = peak which we have broad- 
ened hy a — 0.057r/a. The momentum densities at fc^, = has been 
scaled to the same value, (b) 2D momentum distribution (rife) fol- 
lowing an MI-SF quench for the 2D MI shows a peak atkx = ky = 
which is 3 times higher than the ID peak and qualitatively resembles 
the peak seen in experiments |lj|. 



simple-cubic lattice 112011 which estimates the equilibrium tem- 
perature ~ 2.10J to be below the critical temperature 
Tc = 2.202 J j23|. Since the equilibrium temperature is close 
to the critical temperature, we can use scaling relations around 
the critical point to estimate the momentum distribution and 
the equilibrium temperature below Tc- Using the scaling rela- 
tions for the energy density ll24ll around /3~^ ^ Tc, we obtain 
the equilibrium temperature = Tc— 0.077 J. The momen- 
tum distribution in the vicinity of (3^^ < Tc has a universal 
form at small k 11221 12411 which is given by 



(10) 



Since the critical fluctuation contribution to Uk is only accu- 
rate near A; ^ 0, we have added a ^-independent background 
/o = [1 - w(T)2(l -f CT/27r)] to restore the total spectral 



is the conden- 



weight to unity. Here = B 
sate fraction of the system. The parameters in Eq. [10] have 
been calculated to be: the critical amplitude B = 1.245 , the 
critical exponent (3q w 0.35, the transverse correlation length 

^j,(T) = ^(-) (2^) with ~ 1-65 and v « 0.67 
ll20l - l24ll . To resolve the singularity at A; ~ we have broad- 
ened the momentum distribution by a width cr = 0.057r/a. 
The resulting singular momentum distribution is plotted in 
Fig. [3] as a function of = k, with ky ~ k^ ^ has a 
peak near fc ^ with total weight of 1 — /o ^ 0.3. 



The recent MI-SF quench experiments fP], in the 3D limit, 
can be thought of as a 3D MI-SF quench where the initial 
state is the ground state of the 3D XY Hamiltonian (Eq. |5]) 
with tunneling ^ turned on only along z. This corresponds to 
a lower initial energy density which in turn leads to a reduced 
equilibrium temperature Tc~ (3^^ = . 1 6 J bee ause of energy 
dissipated in the phase-fluctuations along z. This results peak 
near fc ^ in the momentum distribution with total weight 
1 ^ /o ~ 0.6 which is twice as large as the peak obtained 
from quenching from the 3D MI state. 

Conclusion: We have calculated the formation of peaks 
in the momentum distribution of bosonic atoms in an optical 
lattice following a quench from deep in the MI phase to an 
SF phase within the classical XY model approximation. The 
classical results are in semi-quantitative agreement with exact 
iTEBD calculations for the ID MI. Understanding such MI- 
SF quenches is particularly important in the context of recent 
experiments that use such quenches to reach superfluid-like 
states. From our calculations we find that in ID, 2D and 3D 
Mis, the initially uniform momentum distribution develops a 
peak at momentum k = {) following the quench that is qualita- 
tively similar to the one in experimental time-of-flight images. 
However in ID and 2D the resulting momentum distributions 
are non-singular at fc = and represent phases without long- 
range phase order In contrast, we find a true condensate in 
the 3D MI corresponding to the recent experiments! 1| with 
a divergent momentum distribution at fc = 0. One interest- 
ing way to verify our theoretical predictions is to use the re- 
cently developed quantum gas microscope technique [25] to 
look directly at the time evolution of the MI to SF quenches 
in systems of different dimensionalities. 
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